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POTENTIAL THEORY FOR HYPERBOLIC SPDES^ 

By Robert C. Dalang and Eulalia Nualart 
Ecole Polytechnique Federale de Lausanne 

We give general sufficient conditions which imply upper and lower 
bounds for the probability that a multiparameter process hits a given 
set E in terms of a capacity of E related to the process. This extends 
a result of Khoshnevisan and Shi [Ann. Probab. 27 (1999) 1135-1159], 
where estimates for the hitting probabilities of the (A*', d) Brownian 
sheet in terms of the (d — 2N) Newtonian capacity are obtained, and 
readily applies to a wide class of Gaussian processes. Using Malliavin 
calculus and, in particular, a result of Kohatsu-Higa [Probab. Theory 
Related Fields 126 (2003) 421-457], we apply these general results to 
the solution of a system of d nonlinear hyperbolic stochastic partial 
differential equations with two variables. We show that under stan- 
dard hypotheses on the coefficients, the hitting probabilities of this 
solution are bounded above and below by constants times the (d — 4) 
Newtonian capacity. As a consequence, we characterize polar sets for 
this process and prove that the Hausdorff dimension of its range is 
min(d,4) a.s. 

1. Introduction. In this article, we are interested in the following basic 
problem of potential theory for M'^-valued multiparameter stochastic pro- 
cesses: given E C M*^, does this process visit {or hit) E with positive prob- 
ability? Sets that, with probability 1, are not visited are said to be polar 
for the process and otherwise are nonpolar. One objective is to relate these 
hitting probabilities to an analytic expression which is determined by the 
"geometry" of the set, namely the capacity of the set. Another objective is to 
characterize polar sets for the process. In this article, our main goal is to ad- 
dress these questions for non-Gaussian processes that are solutions to a class 
of nonlinear hyperbolic stochastic partial differential equations (SPDEs) in 
the plane (a wide class of Gaussian processes is also considered). 
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There is a large literature about potential theory for multiparameter pro- 
cesses. For multiparameter processes whose components are independent 
single-parameter Markov processes, Fitzsimmons and Salisbury [7] obtained 
upper and lower bounds on hitting probabilities in terms of a notion of en- 
ergy of a set. This type of multiparameter process arises in the study of 
multiple points of single-parameter processes. Song [26] characterized polar 
sets for the n-parameter Ornstein-Uhlenbeck process on a separable Frechet 
Gaussian space as null Cn,2 capacity sets, where the capacity Cn,2 is defined 
in a variational form. Hirsch and Song [8] obtained bounds for the hitting 
probabilities of a class of multiparameter symmetric Markov processes in 
terms of capacity, also in a variational form. 

However, in these references, the class of multiparameter Markov pro- 
cesses does not readily cover certain basic multiparameter processes such as 
the Brownian sheet or the fractional Brownian sheet. In [10], Khoshnevisan 
developed a potential theory for a class of multiparameter Markov processes 
which includes these processes. 

This article is essentially motivated by the work of Khoshnevisan and 
Shi [11], who obtained bounds for hitting probabilities of the Brownian sheet. 
In particular, ii W = {Wt,t G M^) denotes an M'^-valued Brownian sheet, 
they showed that for any compact subset A of M*^ and any < a < 6 < oo, 
there exists a finite positive constant K such that 

Caprf_2jv(^) < P{3 1 G [a, bf :WteA}<K Ca^p^.^^iA), 

where Cap(^_2Ar denotes the capacity with respect to the Newtonian {d — 2N) 
kernel. The proof of the lower bound is essentially based on estimates of the 
first and second moments of functionals of occupation measures. The upper 
bound uses Cairoli's maximal inequality for A^-parameter martingales as a 
key step. 

In this article, we begin by extending their result to a wide class of R"^- 
valued continuous multiparameter processes X = {Xt,t€R^) that are not 
necessarily Gaussian but that have absolutely continuous univariate and bi- 
variate distributions away from the axes. In Section 2, we give sufficient 
conditions on the density of the process that imply upper and lower bounds 
for the hitting probabilities of X in terms of a given capacity related to 
the process (Theorem 2.4). For the lower bound on the hitting probability, 
we require some positivity of a functional of the density of the process (see 
Hypothesis HI) and an upper bound on a functional of the bivariate density 
of the process (Hypothesis H2). For the upper bound on the hitting proba- 
bility, we require that the process be adapted to a commuting filtration (so 
that Cairoli's maximal inequality can be used) and we need a lower bound 
for the conditional density of the increment of the process given the past 
(Hypothesis H3). As a consequence of Theorem 2.4, we obtain an analytic 
criterion for polarity which is given in Corollary 2.5. 
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As a first application of the general result of Theorem 2.4, we consider 
multiparameter Gaussian processes in Section 3. We give sufficient condi- 
tions on the covariance function of a Gaussian process which imply bounds 
for hitting probabilities in terms of the Newtonian capacity (Theorem 3.1). 
This theorem contains many results that exist in the literature and readily 
applies to multiparameter Gaussian processes such as the Brownian sheet, 
a-regular Gaussian fields, the Ornstein-Uhlenbeck sheet and the fractional 
Brownian sheet. For the second and the fourth process, we obtain only a 
lower bound. The upper bound cannot be obtained from Theorem 3.1 since 
such processes are not necessarily adapted to a commuting filtration. 

In Section 4, we apply results of Malliavin calculus and, in particular, the 
very recent result of Kohatsu-Higa [13] to the system of nonlinear hyperbolic 
SPDEs, 

pp. yi d rfiW^ 

(1,1) aMS-g-j^aMS + '-'l-^')' * = 

Xl = xq if tit2 = 0, 1 < i < d, 

where W = {W^ ,j = 1, . . . , d) is a two-parameter d-dimensional Wiener pro- 
cess, the second-order mixed derivative of is the white noise on the 
plane and cj*-, 5* are smooth functions on R'^. It is known [21] that (1.1) has 
a unique continuous solution X = {Xt,t G R^)- In this article, we consider 
this system of equations in the integral form (4.7) as it was studied in [21]. 
In the case 6 = 0, under some regularity and strong ellipticity conditions on 
the matrix a (Conditions PI and P2), we give in Proposition 4.13 an upper 
bound of Gaussian type for the conditional density of an increment of X 
given the past. This uses existing techniques of Malliavin calculus that are 
adapted to the present context (cf. [17] and [18], Chapter 2). We then use 
the result of Kohatsu-Higa [13] to establish a Gaussian-type lower bound for 
the density of the random variable Xt for any t away from the axes (Theo- 
rem 4.14) and use a Gaussian- type lower bound for the conditional density 
of the increment of the process given the past (Theorem 4.16). 

In the last section, we apply the results obtained in Sections 2-4 to the 
solution of system (1.1). In the case 6 = 0, we prove in Theorem 5.1 that 
under Conditions PI and P2 introduced in Section 4.4, the hitting probabil- 
ities of the solution can be bounded above and below in terms of the (d — 4) 
Newtonian capacity. The verification of Hypothesis H2 uses the upper bound 
of Gaussian type obtained in Section 4.4. The main effort in proving Theo- 
rem 5.1 lies in verification of Hypothesis H3, which uses the lower bounds of 
Gaussian type for the density of the solution obtained in Section 4.5. These 
Gaussian-type lower bounds also imply the positivity of the density of the 
solution and so the verification of Hypothesis HI. We treat the case 6^0 
via a change of probability measure (see Corollary 5.3). As a consequence 
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of Corollary 5.3, we prove in Corollaries 5.4 and 5.5 that polar sets for the 
solution to (1.1) are those of (d — 4) Newtonian capacity zero, that the Haus- 
dorff dimension of the range of the solution is min((i, 4) almost surely and its 
stochastic codimension is (d — 4)"'". Finally, we identify d = 3 as the critical 
dimension for the solution to hit points in M"^ (see Corollary 5.6). 

Notice that we obtain the same zero capacity condition for polarity and 
Hausdorff dimension obtained by Dynkin [6], LeGall [16] and Perkins [24] 
for super-Brownian motion. However, there is no direct connection between 
this work and theirs. 

2. General theory. For any s, t G M^, we write s <t when Si < ti for all 
i = 1, . . . , N , where Si denotes the ith coordinate of s, and we write s <t 
when s <t and s ^t. By s At, we mean the point whose ith coordinate is 
Si A for alH = 1, . . . , A^. If s, i e with s<t,we write [s, t] = l\f=i [si, k] 
for an A^-dimensional rectangle and {s,t] =Ylj^i{si,ti]. Finally, we denote 
by II • II the Euclidean norm. 

Let {Q,Q,F) be a complete probability space and let = {!Ft,t EM;!^) be 
a complete, right continuous, commuting filtration, that is, an increasing 
family of sub-cj-fields of Q such that: 

(i) J-Q contains all the null sets of Q. 

(ii) For every i G M^, = n,>t J's- 

(iii) For every s,f G and for all bounded, .Ff-measurable random vari- 
ables y. 



Note that when N = 2, (iii) is hypothesis (F4) of Cairoli and Walsh [1]. For 
N > 2, hypothesis (iii) appears in [10], Chapter 7, Section 2.1. 

Let X = {Xt,t S M^) be a continuous M'^-valued stochastic process defined 
on {Q,Q,¥) and not necessarily adapted to We suppose that for all s,tG 
(0, +oo)^ with t^s, the distribution of the random variable {Xt,Xs) has a 
density that is denoted Pxt,Xs{^^y)- We write pxt{x) for the density of the 
random variable Xt for all t € (0, +oo)^. 

Given a Borel subset E of W^, we denote by V{E) the collection of all 
probability measures on R"' with support in E. 

Following [10], Appendix D, we say that k{-) is a kernel in Mf^ (or a gauge 
function) if k{-) is an even, nonnegative and locally integrable function on 
R"^ which is continuous on R'^ \ {0} and positive in a neighborhood of the 
origin. Basic examples of kernels are the Newtonian f3 kernels A:/3(-) (see 
Section 3.1). 

Given a kernel k{-), for any fi G V{E), we write 



E[Y\J's]=nY\:FsM] 



a.s. 
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and term this quantity the k energy of fj,. The k capacity of E is defined by 

1 



Capfc(S) 



infAtG7'{£;)'^fc(M) 



The following properties of Cap;j(-) are given in [10], Appendix D, in par- 
ticular, Lemma 2.1.2 there. 

Lemma 2.1. Cap^(-) has the following properties: 

(a) Monotonicity. For any two Borel subsets Ei C E2 ofW^, Cap^(£'i) < 
Capfc(i^2). 

(b) Outer regularity on compact sets. For any sequence E, Ei, E2, ■ ■ ■ of 
compact subsets ofR.'^ such that En i E, lim„^oo Cap = Cap^.(£^). 

Let Vo{E) denote the collection of all probability measures on R'^ with 
support in E that are absolutely continuous with respect to Lebesgue mea- 
sure. The absolutely continuous capacity Cap^(ii^) of E with respect to k{-) 
is defined by 

CapO(£;) 



inf;,GT'o(s)'^fc(^)' 

Since Cap^, is not outer regular on compact sets (cf. [10], Appendix D, 
Section 2.2), we define, for all bounded Borel sets E C M.'^, 

Capl^iE) = inf{Cap^(F) :FdE,F bounded and open}. 

Then Cap^(yl) > Cap^'^(^). We now state an additional condition on k 
(which is related to the classical notion of balayage; see [15], Chapter IV) 
that ensures that capacity and absolutely continuous capacity with respect 
to k agree on compact sets. 

Following [10], Appendix D, we say that a kernel A;(-) on M"^ is proper if, 
for all compact sets A C M*^ and fi G ^{A), there exist bounded open sets 
Ai,A2, . . . such that: 

1. An I A. 

2. For all large n > 1, there exist absolutely continuous measures fin with 
support contained in An such that, for all e > 0, there exists A^o such that 
for all n>No: 

(a) fin{An) > 1 - e; 

(b) /iRd k{x - y)iin{dy) < Jj^a k{x - y)fi{dy) for all x £ W^. 



Proposition 2.2 ([10], Appendix D, Theorem 2.3.1). Let k{-) be a 

proper kernel in M.'^. Then, for all compact sets A C M'^, Cap^(j4) = Cap^'^(^). 
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We now introduce the following hypotheses, which ensure a lower bound 
on hitting probabilities for X [see Theorem 2.4(a)]. 

Hypothesis HI. For allO <a<b<oo and M > 0, there exists a finite 
positive constant Ci{a,b, M) such that for almost all x G M"^ with \\x\\ < M, 

I pxAx)dt>Ci. 

Hypothesis H2. There exists a proper kernel k{-) in M.'^ such that for all 
< a < b < oo and M > 0, there exists a finite positive constant C2{a, b, M) 
such that, for almost all x,y €zW^ with \\x\\ < M and \\y\\ < M, 

/ / PXt,Xsix,y)dtds<C2k{x -y). 

J[a,b]'^ J[a,b]^ 

In the case where X is adapted to J^, for s G (0, +00)''^, let Ps{u!, •) be a 
regular version of the conditional distribution of the process {Xt — Xg , t G 
\ [0,s]) given J='s. If for almost ah uj and ah i e \ [0,s], the law 
of Xt — Xg under ^^(0;, •) is absolutely continuous with respect to Lebesgue 
measure on W^, we let Ps,t{^, x) denote the density of Xt — Xs under Ps{u>, •)• 
In this case, there is a null set Ng £ Tg such that for uj £Vl\N E a Borel 
subset of M'^ and s < t, 

G^})= / ps^t{^,x)dx. 

JE 

In particular, ps^ti^^x) is a version of the conditional density of Xt — Xs 
given Ts- The function {uj^t,x) ^ Ps^ti^,x) can be chosen to be measurable. 

Proposition 2.3. Let / : M'^ x ^ M be a nonnegative Borel function, 
let Y be an J^g -measurable random variable and suppose that E[/(Xf — Xg, 
Y)] < 00. Then 

E[fiXt-Xs,Y)\J's]= [ f{x,Y)ps,t{uj,x)dx a.s. 

Proof. If f{x,y) = /i(x)/2(y), then using [5], Theorem 10.2.5, we have 
E[f{Xt - Xs,Y)\J^s] = E[fi{Xt - Xs)f2{Y)\Ts] 
= f2{Y)E[f,{Xt - Xs^ 
= f2{Y) I fi{x)ps,t{u},x)dx 

= / f{x,Y)ps^t{i^,x)dx 



a.s. 
a.s. 
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One easily concludes the proof using a monotone class argument ( [3] , Chap- 
ter I, Theorem 21). □ 

We now introduce a third hypothesis, which leads to an upper bound on 
hitting probabilities for X [see Theorem 2.4(b)]. 

Hypothesis H3. For all < a < b < oo and M > 0, there exists a finite 
positive constant C^{a,h,M) such that for all s G [a,b]^ , a.s., for almost all 

/ Ps,t{uJ,x)dt > C3k{x)l{ii^^Xs\\<M,\\X4<M}{^), 

where k{x) is the same kernel as in Hypothesis H2. 

We are now ready to state the main result of this section. 

Theorem 2.4. (a) Assuming Hypotheses HI and H2, for all < a < 
b < oo and M > 0, there exists a finite positive constant Ki{a,b, M) such 
that for all compact sets ^ C {x G M'^ : ||x|| < M}, 

Ki Capfc(yl) < P{3 1 G [o, bf : Xt G A}. 



(b) If {Xt,t G M^) is adapted to a commuting filtration {J-t,t G M.^), and 
Hypotheses H2 and H3 hold, then for all < a < b < oo and M > 0, there 
exists a finite positive constant K2{a,b, M) such that for all compact sets 
Ac{xG]R'^:||x|| <M}, 

¥{3te[a,b]^ ■.XteA}< K2 Capfc(A). 



Before proving Theorem 2.4, we mention an important consequence. Re- 
call that a Borel set E C M'^ is said to be polar for the process X if 

F{3t e {0, +oo)^ -.XteE} 



0. 



Corollary 2.5. For a process X adapted to a commuting filtration, 
under Hypotheses H1-H3, a compact subset E ofW^ is polar for X if and 
only if CsiY>^{E) = 0. 



Proof. If E is polar for X, then clearly Cap;j(£') = by Theorem 
2.4(a). Conversely, suppose Capfc(£') =0. Write (0, -|-cxd)^ = UmGN[m' "^]^- 
By Theorem 2.4(b), for all m > 1, there is K2 < 00 (depending on m) such 
that 

N 



3t€ 



1 



, m 



m 



-.XteE^^ <K2Capi,iE) = 0. 
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Since this holds for all m, E is polar for X. □ 

Proof of Theorem 2.4. (a) The lower hound. Suppose < a < & < oo 
and < A/ < oo are fixed. Let ^ be a compact subset of {rr G M"^ : ||x|| < M}. 
For e G (0, 1), define = {x : dist(rE, A) < e}, the closed e enlargement 
of A, where dist(x,^) = ||x — proj^x|| and proj^x denotes the orthogonal 
projection of x on A. Fix e G (0, 1) and let / be a probability density on M.'^ 
whose support is contained in A^. We consider the functional Ja b{f) defined 
by 



Ja,b{f)= f{Xt)dt. 
By the Cauchy-Schwarz inequality, 

(2.1) P{3 1 e [a, bf : X, G A} > ^UaM) > 0} > " 
Using Fubini's theorem, we easily deduce from Hypothesis HI that 

(2.2) nJa,b{f)] > Ci 
and from Hypothesis H2 that 

(2.3) E[{Ja,bif)}^]<C2£kif), 

where £k{f) denotes the k energy of the measure /(x) dx. Applying (2.2) and (2.3) 
to (2.1), we obtain 

F{3te[a,bf:XteA}> 



C2£k{f) 

Take the supremum over all /(x) dx G V^A,;) and see that for all e > 0, 
F{3 1 G [a, bf ■.XteA}>^ Capr (^4,). 

By Proposition 2.2, we can replace Capl!^{As) by Cap;j(^£) because k is 
proper. As e — > 0+, | ^, which is compact. By Lemma 2.1(b), Capi^{As) 
converges to Cap^(A) as e — > 0+. Finally, since A is compact and the process 
Xt is continuous, 

fl {3 1 G [a, 6]^ : G = {3 1 G [a, 6]^ : G ^}. 

£>0 



We conclude that 



P{3t G [a,hf:Xt G A} > ^Capfe(A). 



C2 

This concludes the proof of (a) of Theorem 2.4. 



POTENTIAL THEORY FOR HYPERBOLIC SPDES 9 

(b) The upper bound. Suppose < a < 6 < oo and < -A/ < oo are fixed. 
Let ^ be a compact subset of {x G M'^: ||x|| < M}. Let /iM'^ ^ M+ be a 
measurable function such that E[{ Jf,^2fe-a(/)}^] < oo. We define the following 
square integrable multiparameter martingale: 

Mt(/)=E[Jb,2fe-a(/)|^t], 

Since J^t is a commuting filtration, we can apply Cairoli's maximal inequality 
(see [10], Chapter 7, Theorem 2.3.2), to get 



E 



sup {Mtif)Y 

iG[a,b]^nB^ 



< 4^ sup E[{Mi(/)}2] < 4^E[{ J,,26-a(/)}' 



where B denotes the set of dyadic rationals. Suppose that / is a density 
function on supported on {x S M*^: < A/}. By Hypothesis H2 and 
(2.3), we get 



(2.4) E 



sup {Mtif)y 

ie[a,6]^nID)JV 



<4^C2£k{f)- 



For e G (0,1), define = {x G M'^ : dist(2;. A) < e} , the open e enlargement 
of A. Suppose e is small enough so that AeC{xe W^:\\x\\ < M}. We can 
assume that 

(2.5) F{3t£[a,b]^ :Xt€Ae}>0 for ah e > 0. 

Indeed, if there exists an e > such that this probability is equal to zero, 
then the upper bound is trivial since 

¥{3t G [a, b]^ ■.Xt£A}<F{3tG [a, b]^ : Xt G As} for all e > 0. 

Assuming (2.5), we claim that there exists a random variable T^, taking 
values in ([a, b]^ D D^) U {+00} , such that 

(2.6) {r^<oo} ^ {3te[a,bf DB^ -.XteAs} 

and Xt^ G A^ on {T^ < 00}. Indeed, order the set [a, b]^ Pi = {qi,q2, ■ ■ ■} 
and define = qin{{k:X eAeji where inf is defined to be +00 and in this 
case = +00. Note that assumption (2.5), the fact that A^ is open and the 
continuity of 1 1— > imply that ¥{T^ < 00} > 0. In particular, it is possible 
to condition on the event {T^ < 00}. 
For any Borel set E CW^, define 

fie{E) = ¥{Xt^ G E\T^ < 00}. 

Clearly G V{As). Moreover, is absolutely continuous because every Xf 
is: Let fe{x) be the Radon-Nikodym derivative of fie with respect to A. Then 
fe is supported on A^. Applying (2.4) to fe, we have that for all e G (0, 1), 



E 



sup {MtUe)Y 



<4^C2£k{fe 



10 R. C. DALANG AND E. NUALART 

We claim that for all e G (0, 1) and all s G [a, 6]^, 

(2.7) MsUe) > C3l{||X.||<Af} / fe{x + Xs)k{x) dx 



a.s. 



Indeed, fix s E [a,b]^ . Using the fact that Xg is ^s-measurable, Proposi- 
tion 2.3 and Hypothesis H3, we have that for almost all w S 



Msife 



f fe{Xt-Xs+Xs)dtJ 

/ fe{x + Xs)( ps,t{uJ,x)dt] dx 

JR'i \J\b.2b-a\ J 



a.s. 



> C'3l{||X,||<A/} / je{x^Xs)k(x)dx 



a.s. 



The last inequality follows since is nonnegative and supported in . This 
completes the proof of (2.7). 

Since (2.7) holds for all s G [a,&]^ and since G [a, 6]^, we can replace 
s by in (2.7). Note that < M} holds on {T^ < cx)}. Therefore, 

(2.8) sup Msife) > Csl^T^^oo} f fe{x + XT^)k{x)dx a.s. 

Square both sides of the last inequality, take expectations and apply (2.4) 
to the left-hand side, to obtain 



4^C2£k{fs)>ClF{T'<oo}E \^j^Je{x + XTe)k{x)dxj <oo 

= ClnT'<^}l (I h{x + y)k{x)dx\%,{y)dy. 
Using Jensen's inequality, we get 

^^C2£kife) > ClF{T' < oo} f / f k{x- y)fe{x)fe{y) dx dy 

= ClF{T'<oc}{£kife)f. 



If £k{fe) were finite, this would imply 



(2.9) 



»{r^ <oo}< 



CiSkifeV 

but we do not know a priori that fs has finite energy. For that reason, we 
use a truncation argument. 

For all g > and all e £ (0, 1), define 

f^{x)=Mx)l[0,,]{fe{x)), xeR". 
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Since fs is supported on A^, so is f^. Moreover, the latter is a subprobability 
density function that is bounded above by fe and q. Therefore, since k is 
locally integrable in W^, <?fc(/|) < oo. Apply to f§ exactly the same argument 
that led to (2.8) to see that 

sup M,(/|)>C3l{T.<oo} / f^{x + XT^)k{x)dx a.s. 

Square both sides of the inequality, take expectations and use Jensen's in- 
equality to get 



E 



sup {M,(/|)}2 

3e[a,fe]^nD^ 



>C|P{r^<oo}f / / k{x-y)f^{x)fe{y)dxdy)\ 

By (2.4), the left-hand side is bounded above by 4^C2<?fc(/|). The right- 
hand side is clearly bounded below by C^F{T'^ < oo}(<5fe(/|))^. Hence, we 
obtain 

F{T' <oo}< ^ 



I)- 



Finally, since k{x) is nonnegative, lim^i+oo '?A:(/|) =£k{fe)i we can let q1 
+00 in the above inequality and use the monotone convergence theorem to 
obtain (2.9). 

Now, since t*-^ Xt is continuous and is open, using (2.6) and (2.9), we 
obtain 

F{3t£[a,bf ■.XteAe}< ' 



CiSkifJ-e)' 

Recall that fis G T'(^e), so for all e > 0, 

P{3 1 G [a, bf : G ^} < P{3 1 G [a, : Xt G A^} < Capfc(A,). 

Finally, since A^r is compact, Lemma 2.1(b) implies that Cap;i.(^£) converges 
as e ^ 0+ to Cap;,(j4). This concludes the proof of Theorem 2.4(b). □ 



3. Multiparameter Gaussian processes. In this section, we focus on Gaus- 
sian processes and reformulate the Hypotheses H1-H3 as conditions on the 
covariance of the process so as to relate bounds on hitting probabilities di- 
rectly to properties of the covariance. 



12 



R. C. DALANG AND E. NUALART 



(3.2) 



1 



<C3||t-s|r if\\t-s\\<6, 



3.1. Relating Newtonian capacity and covariance. Let (fi, ^, P) be a com- 
plete probability space. Let X = {Xt,t € M^) be a continuous R'^-valued 
centered Gaussian process with independent coordinate processes (Xl). For 
all t G (0,00)''^, we denote by pxti^) the density function of the centered 
Gaussian random variable Xt on M*^. For all s, t G (0, 00)^, we write a{s, t) = 
K[XlXf], which does not depend on i, cr^(t) = (j{t,t) and p{s,t) = a{s,t)/{a(t)a{s)). 

Given a £ (0, 1) and 7 > a, we introduce the following hypotheses. 

Hypothesis Al. For all < a <b < 00, there exist positive finite con- 
stants 6, e and Ci, . . . , C5 such that for all s,t£ [a, h]^ , 

(3.1) Ci<a'^{t)<C2, 

<^{s,t) 
a^s) 

(3.3) C^t-sf" <l-p^{s,t)<C5\\t-sf'' if\\t-s\\<S, 

(3.4) \p{s,t)\<l-e if\\t-s\\>5. 

Hypothesis A2. For all < a < b < 00 and M > 0, there exist finite 
positive constants Cq, C-j and Cg such that for all s,t£ [a, 6]^ with s <t, 

(3.5) l^\\Xs\\<M}\nXt - XslJ's]] < Ce\\t - s\\\ 

(3.6) CjWt - " < E[{Xt - Xs) - E[Xt - Xs\J^s]f < CsWt - sf". 

Theorem 3.1. (a) Assume there are a € (0,1) and 7 > a for which 
Hypothesis Al holds and N/a > 2. Then for all < a <b < 00 and M > 0, 
there exists a finite positive constant Ki{a,b, M) , such that for all compact 
sets Ac{x£R'^: \\x\\ < M}, 

K, Cap,_(^/„) {A)<¥{3te [a, bf:XtGA}, 

where, for /3 > 0, Cap^(-) denotes the capacity with respect to the Newtonian 
j3 kernel kp{-), where 




ifO<(3<d, 
if (3 = 0, 



an 



kp{x) = 

dfor (3<Q, Cap^(-) = 1. 

(b) Suppose {Xt,t G M^) is adapted to a commuting filtration {J-t,t G 
M^). Assume there are a G (0, 1) and 7 > a for which Hypotheses Al and 
A2 hold and N/a > 2. Then for all < a <b < 00 and M > 0, there exists a 
finite positive constant K2{a,b, M) such that for all compact sets Ac {x G 
R'^:\\x\\ <M}, 

P{3 1 G [a, bf ■.XteA}<K2 Cap,_(jv/„) (A). 
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Remark 3.2. (a) For < (3 < d, the functions kjs are not only posi- 
tive but even positive-definite: this follows from [27], Chapter V, Section 
1, Lemma 2(b). For /3 = 0, the kernel is not nonnegative on R'^, but 
koix — y) > for x,y ^ A when Ac {x GM.'^: \\x\\ < M}, and this is suffi- 
cient for the results of Section 2 to hold. 

(b) Note that for any < /3 < d — 2, A;^(-) is a proper kernel. Indeed, given 
a compact set Ac {x : \\x\\ < M} and /j, € 'P{A), for any n > 1, let An = 
{x G : dist(x, A) <1 /n} be the open enlargement of A. Let {Bt, t>0) be 
a standard Brownian motion in M'^ and let pt be the density of Bt. For n > 1, 
let fin,t be the measure whose density function is the restriction to the set An 
of Pt * where * denotes the convolution product. Set /(x) = (// * kp){x). 
Since AA:/3(x) < for < /? < d - 2, E^{f{Bt)) < f{x) for ah t > and 
X E M*^, or equivalently. 



kp{x-y){pt* fJ'){y)dy< / kp{x-y)fi{dy). 

JR'' 

For small t, the fin,t ^ire nearly probability measures, and so (a) and (b) of 
the definition of a proper kernel hold. 

(c) The condition N/a > 2 ensures that kd~[N/a) is a proper kernel. This 
is only a restriction when = 1 . 

(d) If a = 1/2 and d < 2N, the choice Cap^_2jv(^) = 1 is natural, since 
in this case, the Brownian sheet hits points in R'^ (cf. [23]). 

Before proving Theorem 3.1, we mention an important consequence. Given 
s > and a Borel subset E of M"^, let 

{oo oo ^ 

Y,i2r,y:Ec\jB{xi,ri),snpri<e\, 
i=i i=i « J 

where B{x, r) denotes the closed ball of radius r > centered at rr G M*^. The 
TCs is called the d- dimensional Hausdorff measure. Moreover, we associate 
to the set E a number dim7^(£') as follows: 

dim^^ (E) = sup{s > : (^) = oo} = inf {s > O-.HsiE) = 0}. 

This is the Hausdorff dimension of E. Following [10], the stochastic codi- 
mension of a random set E in W^, denoted codim(^^), if it exists, is the real 
number /? E [0,d] such that for all compact sets A C W^, 



'{EnA^0} 



> 0, whenever dim>^(yl) > /3, 
= 0, whenever dim'j-c{A) < j3. 



The following result gives a relationship between Hausdorff dimension and 
stochastic codimension. 
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Theorem 3.3 ([10], Theorem 4.7.1, Chapter 11). Given a random set 
E in whose codimension (5 is strictly between and d, 

dim-}-i(E) + codim(£') = d a.s. 

Theorems 3.1 and 3.3 imply the fohowing result. 

Corollary 3.4. Under the hypotheses of Theorem 3.1(b), 
codim{X((0,+oo)^)} = {d- (iV/a))+ 

and if d> N/a, 

dimn{X{{0,+oo)^)} = N/a a.s. 

Proof. By Frostman's theorem (see [10], Appendix C, Theorem 2.2.1), 
the capacitarian and Hausdorff dimensions agree on compact sets. Therefore, 
the desired result follows from Theorems 3.1 and 3.3. □ 

Proof of Theorem 3.1. (a) Suppose < a < 6 < oo and M > are 
fixed. We assume that Hypothesis Al holds for a G (0, 1) and 7 > a fixed, 
and show that Hypotheses HI and H2 of Theorem 2.4 hold. 

Verification of Hypothesis HI. Fix x G M'^ such that \\x\\ < M. 
Inequality (3.1) implies that 

/ pxAx)dt = {2n)-''/' [ a-'{t)e^p(-j%^dt 

Jla,b]^ J[a,b]N \ 2a^{t)J 

>(2vr)-^/^(C,)-^/^(6-a)^exp(-^), 

which proves Hypothesis HI. 

Verification of Hypothesis H2. By (3.3) and (3.4), for ah s,tG 
[a,b]^ with t / s, {Xt,Xs) has a (Gaussian) density PXt,Xs{x,y). The latter 
can be written as 

PXt,Xs{x,y) =PXt\Xs=y{x)VXs{y), 

where Pxt\Xs=y denotes the conditional density function of the random vari- 
able Xt given Xs = y. 

Note that the conditional distribution of XI given X] = is Normal with 
mean m{s,t)y'^, where m(s,t) = {a{s,t))/{a'^{s)), and variance r^(s,t) = 
o-2(i)(l - p2(s,t)). Observe that r2(s,t) > by (3.3) and (3.4), that (3.2) 
is a condition on the conditional mean m(s,t) and that (3.3) is a condition 
on the conditional variance. 
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Fix x,y£ such that ||x|| < M and ||y|| < M. By (3.4), pxt,x,(-,-) is 
bounded by some constant C" > when ||t — s|| > 6. Therefore, 

/ / PXt,Xs{x,y)dtds<C + pxt,Xs{x,y) dtds, 

J[a,b]N J[a,b]N J JD{5) 

where D{6) = {{s,t) G [a,b]^ x [a, 6]^ : ||t - s|| < ,5} and C = C'{b- a)^^. 
The integral on the right-hand side can be written 

\\x — m{s, t)y\\' 



(27r)"'^// r"'^(s,t)exp 

(3.7) 



D{S) ' V 2r2(s,t) 



X cr "^(g) exp( — ]dtds. 

* 2(7^ (s) ' 



By the triangle inequality and the identity {u — v)'^ > — f^, 
exp 



(3.8) 



/ \\x -m{s,t)y\\^ \ 



<exp(-M%xprll^ll^l^— 



4r2(s,t); 2r2(s,t) 
By (3.2), there exists a constant C3 such that 

(3.9) |l-m(s,t)p<C3p-sfT. 
By (3.1) and (3.3), 

(3.10) Ki\\t-sf" <T\s,t)<K2\\t-sf''. 

We now apply (3.8)-(3.10) to (3.7). Because 7 > a, this yields 

PXt,Xs{x,y)dt ds 



[a,b]N J[a,b]N 

<C+(2^)-^(iri)-^/2(^^)-d/2 f r dtds\\t-s\\-'''^ 

J Jd(S) 



<C + K3ll \\t - 5||-"'^exp( - ^ J^,^ ^ ) dtds. 



2 

We now fix t and use the change of variables u = t — s to see that this 
expression is less than or equal to 
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where B{6) = {u :\\u\\ < 6}. Finally, use the change of variables u ■ 
\\x — y\\^^"z{4K2)~^^^'^"^ to see that this is less than or equal to 



C + K4\\x-y\ 



-d+iN/a) 



X / \\z\\-'''^ exp{-l /\\zf")dz. 



I B{{4K2y/'-^°''>S/\\x-y\\^/") 

We now state a real variable technical lemma which is crucial for our 
estimates. The proof of this lemma is left to the reader. 

Lemma 3.5. Define ipa^pir) = /^^.^^ ||z||~^e~^/ll^ll^" dz, /or a// r > 0. Then 
for any rg > and a G (0, 1), there exist finite constants 01,02,03, C4 > such 
that for all r > rg, 

ci<^Pa,/3{r) <02, if(3>N, 
C3ln(r/ro) < ^a,p{r) <Ci\n{r), if [3 = N. 

Continuing the verification of Hypothesis H2, apply Lemma 3.5 with 
/3 = adto (3.11) and use the fact that C < C(2M)^-W") be- 
cause < M and ||y|| < M, to conclude the verification of Hypothesis H2 
for d > N/a. When d = N/a, choose ro > such that (4K2)^/^^"^'5/(2Af )^/" > 
ro and apply Lemma 3.5 to (3.11) to obtain 

-(4^)V(2«)(^^ 



/ / PXt,Xsix,y)dtds <C + 0iK4ln[— Trrn;-]- 



Note that for all x,y G M'^ with ||3;|| < M and ||y|| < M, y|| < 2M. Then 
we can check that if ||a; — ?/|| < 2M, there exists a finite constant C" > 1 such 
that 



y V ~ y\ 

On the other hand, note that C < C(ln(3/2))"Mn(3M/||x - y||), and the 
verification of Hypothesis H2 for d = N/a is completed. When d < N/a, the 
expression in (3.11) is bounded, so Hypothesis H2 holds with k{x) = 1. 

This completes the proof of Theorem 3.1(a). 

(b) We now assume that Hypotheses Al and A2 hold for a S (0, 1) fixed 
and 7 > a, and show that Hypothesis H3 of Theorem 2.4 holds. We also 
assume that d > N/a, since otherwise, the statement is trivial. 

Verification of Hypothesis H3. For all s,te (0,oo) with s <t, and 
for almost all x E M*^ and u; S let ps^t{^,x) denote the conditional density 
of the random variable Xt — Xg given Tg- By (3.4), the latter exists and. 
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for almost all u £ is a Gaussian density with conditional mean ^{s,t) = 
E[Xt - Xs\Ts] and deterministic variance /?^(s,t) = E[(Xj - Xs) - E[Xt - 
Xs\Ts]? (see [4], Chapter II, Section 3). 

It suffices to check the inequality in Hypothesis 113 when ||x|| < 2M, since 
the indicator on the right-hand side vanishes for ||x|| > 2M. Fix s E [a, 6]^. 
Then 

2/32(s,t) 

By the triangle inequahty, 

(3.13) exp —J- — - — >exp -— ^- — — exp' 



f3 12) -^[^'^f--"]^ „ 

={2.)-"'i p-\s^t)..j j^-j:i'^f ']dt. 



2/3^s,t) ;-^^^V 2(3^s,t)J ' \ 2l3^s,t) J' 

We now apply (3.5), (3.6) and (3.13) to (3.12), and we see that this 
expression is greater than or equal to 



Using the fact that 7 > a and the change of variables t — s= \\x\\^^°'z{C'i)~^^^'^'^^ , 
we get that this is greater than or equal to 

(3.14) i^5||x|r'^+W") / ||z|r"'^exp(-l/||zf")dz. 

J Bi(C7y/(^'^) {b-a)/\\x\\^/'^) 

It now suffices to choose tq > such that {C'r)^^^^'^\b - a)/{2My/'^ > 
3^/"ro and apply Lemma 3.5 to (3.14). This concludes the verification of 
Hypothesis H3 when d > N/a. Finally, for d = N/a, we apply Lemma 3.5 
to (3.14), to obtain 

f / N , ^ /6i/"Mi/"\ 1 , /3M\ 

/ Ps,t{uJ,x)dt>C3K5ln[ > -csAsln — — , 

which concludes the proof of Hypothesis H3 for d > N/a. 
This completes the proof of Theorem 3.1(b). □ 
3.2. Examples. 

3.2.1. The Brownian sheet. Suppose that W = {Wt = (W/, . . . , Wf),t G 
M.^) is a d-dimensional A^-parameter Brownian sheet, that is, the coordinate 
processes of W are Gaussian, with zero mean and covariances 

N 

E[W^W^] = Y[{siA ti)5kj for all s, t G M^, 1 < k, j < d, 

i=l 



X\ 



2 \ ^ ci\\t-sp^ 
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and defined on the canonical probability space (O, ^,P), where 0, is the space 
of all continuous functions oj : — > M.'^ vanishing on the axes, P is the law 
of W and Q is the completion of the Borel c-field of 0, with respect to P. We 
also consider the increasing family of <t- fields J- = {J-t,t G M^), such that 
for any t G M^, is generated by the random variables (Ws,s < t), and 
the null sets of Q. The latter is a complete, right continuous, commuting 
filtration (see [10], Chapter 7, Theorem 2.4.1). 

Theorem 3.6 ([11], Theorem 1.1). The {N,d) Brownian sheet satisfies 
Hypotheses Al and A2 with a = 1/2 and 7 = 1, and therefore the conclusions 
of Theorem 3.1. 

Proof. Fix < a < b < oo. Inequality (3.1) is trivial since E[W(^] = 
nili ti, for ah t e M^. Moreover, since the Brownian sheet has independent 
increments, (3.5) holds also trivially. 

We claim that for all s,t G [a, 6]^ with s <t, there exist two constants 
Ci and C2 such that 

(3.15) Ci||t - s\\ < E[{Wt - Wsf] < C2\\t - s\\. 

Suppose s<t and note that E[{Wt - Ws)'^] = a'^{t) - a'^{s) = Uf=iti - 
nili Si. Let f{t) = nili ti and let (j) = {<j)i, . . . , c^n) : [0, 1] be an affine 

function such that (p{0) = s and = t. The partial derivatives of / are 
bounded above and below on [a, 6]^, so 

fit) - f{s) = fm)) - fim) = / E 7^(0N)'A.(n) du 

JO UTi 

N 

<cY,m'^)-(t>m)<c'\\t-si 

i=l 

which proves the upper bound in (3.15). The lower bound is obtained by 
proceeding along the same lines. 

We now prove (3.2) and (3.3). If s,t £ [a,h]^ , using (3.1) and (3.15), 

cr2(s) -cr(s,t) _ (s) -0-2(5 At) 



a^{s) 



a^{s) 



<C\\s-{sM)\\<C\\t 



which proves (3.2). To prove (3.3), if s,t G [a,^] , using (3.1) and (3.15), 

a'^{t)a'^{s) — {a{s,t)Y a'^(sAt) ,9,, , , 

Jit)ail ^ 2^)^ - .^(. A t) + a\s) - a^s A t)) 

>C{\\t-{sAt)\\ + \\s-{sAt)\\) 
>C\\t-s\\, 
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which concludes the proof of (3.3). 

Since (3.6) fohows from (3.15), it remains to prove (3.4). Because s,ti— > 
p{s,t) is continuous on [a,b]^ , it suffices to check that p{s,t) < 1 for any 
s,t £ [a, ^^l'^ with s ^t. This is clear since 

p(.,t)=n^<i- 



This proves Theorem 3.6. □ 

3.2.2. a-regular Gaussian fields. For a£ (0,1), an A^-parameter, M*^- 
valued process X = {Xt,t G W^) is said to be a-regular (see [10], Chapter 
11, Section 5.2) if X is a centered, stationary Gaussian process with i.i.d. 
coordinate processes whose covariance function R satisfies the following: 

Assumption R1. If s / 0, then |-R(s)| < 1, and there exist finite positive 
constants ci < C2 and 6 > 0, such that, for all s G with ||s|| < 6, 

ci||sf" < 1-R{s) < C2||sf°. 

Such Gaussian fields were studied in [10], Chapter 11, Section 5.2. In 
the latter, the stochastic codimension and the Hausdorff dimension of the 
range of these processes is obtained (see [10], Chapter 11, Theorem 5.2.1 
and Corollary 5.2.1). 

The following result proves that the a-regular Gaussian fields satisfy the 
lower bound of Theorem 3.1(a). The upper bound cannot be obtained from 
Theorem 3.1(b) since such processes are not necessarily adapted to a com- 
muting filtration. An upper bound involving Hausdorff measure is given 
in [10], Chapter 11, Example 5.3.3. 

Theorem 3.7. Fix a G (0,1) with N/a > 2. Let X = {Xt,t G M^) be 
a-regular. Fix < a < b < oo and < M < oo. Then there exists a finite 
positive constant K such that for all compact sets A C {x G M"' : ||x|| < M}, 

K Caprf„(^/,) {A) < P{3 1 G [a, bf : Xt G A}. 

Proof. Under Assumption Rl, R{0) = 1, so cr^(t) = 1 and p{s,t) = 
R{t — s) = a{s,t). In particular, (3.1) with 7 = 2a and (3.2) hold. Condi- 
tion (3.3) holds (for a sufficiently small (5 > 0) because 

1 - p\s,t) = (1 + p{s,t)){l - pis,t)) > Cl(l - C2||i - S||2°)||i - 
>Ci(l-C2(^2")||t-s||2°. 

Since R{-) is a bounded covariance function, it is nonnegative definite. By 
Bochner's theorem ([25], Theorem 6.5.64), R(-) is the Fourier transform 
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of a nonnegative measure, which is a probabihty measure since R{0) = 1. 
Therefore, R{-) is in fact continuous, so by Assumption Rl, there is e > 
such that \R{t-s)\ < 1-e for > 6 with s,t€ [a,b]^ . Therefore, (3.4) 

holds. This proves that Hypothesis Al holds with 7 = 2a, so the conclusion 
follows from Theorem 3.1(a). □ 

3.2.3. The{N,d) Ornstein-Uhlenbeck sheet. The (A^, d) Ornstein-Uhlenbeck 
sheet U=iUt,t£ M^) is defined as 

where \t\ = \ti\-\ \-\tN\, e* = (e*i , . . . , e*^) and W = {Wt,t gR^) is an M"'- 

valued Brownian sheet. Therefore, U is an A^-parameter centered stationary 
Gaussian process on R"' with covariance function given by 



+ ■ 



Consider its natural and completed filtration, denoted J-jr , that is, 

.Ff = a{Us,s <t}= a{Ws, s < e*} = jj, t G M^, 

where J-]^ denotes the natural filtration of the Brownian sheet and is there- 
fore a commuting filtration. 

Theorem 3.8. The {N,d) Ornstein-Uhlenbeck sheet satisfies Hypothe- 
ses Al and A2 with a = 1/2 and 7 = 1, and therefore the estimates of 
Theorem 3.1. 

Proof. Since the (A^, d) Ornstein-Uhlenbeck sheet is an a-regular Gaus- 
sian field with a = 1/2, the lower bound follows from Theorem 3.7 and it 
remains to verify Hypothesis A2. Using the estimate 1 — < x for all x > 
and using (3.15), we have, for all s,t G [a, with s <t, 

I g I" jy^ jj I jrU j I 

= |E[e-l*l/2(T^,* - Wes) + (e-l*'/2 - e-^'\/^)WeA:F^]\ 
(3.16) =|(e-l*l/2-e-l^l/2)t^,.| 
<\Us\i\t\-\s\)/2 
< CiM||t- s|| 

on {|C/s| < M}, which proves (3.5). To prove (3.6), use (3.16) to see that the 
expectation in (3.6) is equal to 

E[e-l*l/2t^,. - e-l^l/Ve. - (e-l*l/2 - e-l^l/2)VFe»]' = E[e-\'\/\W,. - W^s)]', 

so by (3.15), there exist two constants C2 and C3 such that for all s, t G [a, b]'^ 
with s <t, 

C2\\e^ - e'W < e-l*lE[(Wet - We^f] < Cglje* - e"||. 
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Finally, using the inequalities cx <1 — e ^ < x, for x £ [0,b], we obtain 

C4t - s\\ < e-l*lE[(Wet - Wesf] < Csllt - s\\, 
which proves (3.6). The upper bound follows then from Theorem 3.1(b). □ 

3.2.4. The fractional Brownian sheet. The fractional Brownian sheet with 
Hurst parameter H € (0, 1), = {X^ ,t E M^), is a d-dimensional centered 
Gaussian process with independent coordinate processes, each with covari- 
ance function given by 

N 

E[{Xf^y{X^y] = n ^7:{tr + - \U - s^r\ G M^, 1 < J < d. 



where c is a positive finite constant. Note that if H = ^, one obtains the 
standard Brownian sheet. 

As in Theorem 3.7, we obtain only a lower capacity estimate for the 
fractional Brownian sheet as a consequence of Theorem 3.1(a), since this 
process is not necessarily adapted to a commuting filtration. 

Theorem 3.9. Fix H e (0,1), < a < b < oo and < M < co. Then 
there exists a finite positive constant K such that for all compact sets A C 
{xGM^:||x|| <M}, 

K C^Vd-iN/H) {A) < P{3 1 G [a, hf : X^ G A}. 

Proof. By Theorem 3.1(a), it suffices to prove that Hypothesis Al 
holds with a = and 7 = 2HAI. Fix < a < 6 < oo. Since ^^(t) = cJlili 
for ah t G M^, (3.1) holds trivially. On the other hand, for all s,t€ [a,b]^ , 
we have 



(3.17) 



a{s,t) 



a^s) 



nili - uliiisF + - \h - ^.P^)/2) 



To prove (3.2) for s fixed, set /(n) = nr=i |(s?^ + uf^ - W - Si\^^) and 
let = (01, . . . , (/>7v) : [0, 1] — > be an affine function such that (^(0) = t 
and (/)(1) = s. Because / is differentiable in the orthant centered at s that 
contains and in this orthant, for ||t — s|| sufficiently small and u G Imi;^>, 
|^(n)| < C(\ui — Sip^"^ V 1), the numerator on the right-hand side of (3.17) 
can be bounded by 

1 TV 

|/(,^(1)) - /(0(O))| < C \ E(|</).(r) - V \)4>i{r) dr = C\t - s\''"''\ 

■^0 i=l 

This concludes the proof of (3.2). 
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To prove (3.3), set /(u) = Uf=i9{ui), where g{ui) = ^^-^(1 + u]^ - 
\ui — Ip^)- An elementary calculation shows that for \ui — 1| sufficiently 
small, there are constants C2 > ci > such that ci\ui — < \g'{ui)\ < 

C2\ui — Therefore, there are constants C2 > ci > such that for 

||n — (1, . . . , 1)11 sufficiently small, 

Of , . 

-w 



(3.18) ciliii - Ip^"^ < 



dui 



<C2\ui-l\'^^-\ 



Let ip = ■ • ■ ) 4'n) '■ [0, 1] — > be an affine function such that 
(1, . . . , 1) and (/>(0) = (ti/si, . . . ,tj\[/sN)- Observe that 



1 ^ 



1 - p{s, t) = mi)) - nm) = i E ^('^("))'^^(^) 



i=i 

By (3.18), for s,t G [a, 6]^ with \\s — 1\\ sufficiently small, this expression is 
bounded above and below by constants times 

TV 1 N 
E / mr)-ir-'Mr)dr = Y^\U-s,r/s]^. 
1=1-^0 i=l 

Because 1 — p'^{s,t) = (1 — p(s,i))(l + p{s,t)) and the second factor is less 
than or equal to 2 and bounded away from 0, it follows that (3-3) holds. 

It remains to prove (3.4). For this, it suffices to check that for any Si,ti € 
[a, b] with ti — Si> 0, 



Move the square root to the left-hand side and isolate the perfect square 
to see that this is equivalent to tf — s^^ < \ti — Si\^ , which holds since 
< H < 1. The lower bound follows then from Theorem 3.1(a). □ 

4. Gaussian-type bounds for densities of solutions of hyperbolic SPDEs. 

In this section, we use techniques of Malliavin calculus to establish Gaussian- 
type bounds on the density of the solution to a nonlinear hyperbolic SPDE. 
These bounds are such that they can be used to verify Hypotheses H1-H3, 
as we see in Section 5. 



4.1. Elements of Malliavin calculus. In this section, we recall, follow- 
ing [18], some elements of Malhavin calculus. Let W = {Wt = (W/, . . . , Wf), 
t S R^) be an M'^-valued two-parameter Brownian sheet defined on its canon- 
ical probability space {Q,Q,F) and let = {J-t,t G M^) be its natural ffitra- 
tion (see Section 3.2.1). 

Let H be the Hilbert space H = L'^{Rl,R'^). For any h e H, we set 
W{h) = E1=iIr2 hj{z)dWi. The Gaussian subspace H = {W{h),he H] 
of L^(ri,^,P) is isomorphic to H. 
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Let S denote the class of smooth random variables F = f(W{hi), . . . , W{hn)), 
where hi,. . . ,hn are in H, n>l, and / belongs to Cp'(R"), the set of func- 
tions / such that / and all its partial derivatives have at most polynomial 
growth. 

Given F in S, its derivative is the d-dimensional stochastic process DF = 
{DtF = {d'^^F, DfF),t G ^l) given by 

n 

DtF = jr^iWihi), W{K))hiit). 

i=l 

More generally, the fcth-order derivative of F is obtained by iterating the 
derivative operator k times: if is a smooth random variable and k is an 
integer, set D^^ j-^F = Dt^ ■ ■ ■ Dt^.F. Then for every p>l and any natural 
number k, we denote by n^'P the closure of S with respect to the seminorm 
II • \\k,p defined by 

\\F\\i,=n\Fn+j2n\\D^Fr^^,], 

where H^^ is the product space L'^ {{R\y ,R'^) , and 

\\D^F\\^.. = ( y: ••• 
\ki,...,kj=i-^^+ •^^+ 

We set B- = np>infc>iB^'P. 

Similarly, for any separable Hilbert space V, we can define the analogous 
spaces D'^'P(y) and D°°(y) of ^-valued random variables, and the related 
II ■ llfc.p.y seminorms (the related smooth functionals being of the form F = 
J2]=i FjVj, where Fj e S and Vj gV). 

We denote by 6 the adjoint of the operator D, which is an unbounded 
operator on L?'{Q.,H) taking values in L^(0) ([18], Definition 1.3.1). In par- 
ticular, if n belongs to DomJ, then 5{u) is the element of Lp'iQ) characterized 
by the duality relationship 

(4.1) E{F5{u)) = E [y1^ J ^ d[^^f4 dt^ for any F G D^'^. 

If n G L^(]R^ X iljM'^) is an adapted process, then (see [18], Proposition 1.3.4) 
u belongs to Dom5 and 5{u) coincides with the Ito integral 

d „ 

5{u)=y / 4dwi. 

We use the following estimate of the || • \\k^p norm of 5{u). 
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Proposition 4.1 ([18], Proposition 3.2.1, and [19], (1.11) and page 131). 
The adjoint 6 is a continuous operator from n'''^^'P{H) to D^''^ for all p> 1, 
k>Q. Hence, for all u £ D^+i'P(i7) , 

(4-2) ll^(li)l|fc,p < Ck,p\\u\\k+l,p,H 

for some constant c^^p > 0. 

Our first application of Malliavin calculus to the study of probability laws 
is the following global criterion for smoothness of densities. 

Theorem 4.2 ([18], Theorem 2.1.2 and Corollary 2.1.2). Let F = {F^, 
. . . , F'^) be a random vector satisfying the following two conditions: 

(i) For all i = l,...,d, F*eD°°. 

(ii) The Malliavin matrix of F defined by = {{DF^ , DF^) H)i<ij<d is 
invertible a.s. 

Then the probability law of F is absolutely continuous with respect to Lebesgue 
measure. Moreover, assuming (i) and (ii), and 

(iii) {det -fp)'^ e LP for all p>l, 

the probability density function of F is infinitely differentiable. 

A random vector F that satisfies conditions (i)-(iii) of Theorem 4.2 is 
said to be nondegenerate. For a nondegenerate random vector, the following 
integration by parts formula plays a key role. 

Proposition 4.3 ([19], Proposition 3.2.1). Let F = {F^, . . . , F'^) e (0°°)*^ 
be a nondegenerate random vector, let G G and let g E C^{W^). Fix 
k > 1. Then for any multiindex a = (ai, . . . , a^) G {1, . . . , d}^ , there exists 
an element Ha,{F,G) G Ji°° such that 

(4.3) n{dag){F)G]=n9{F)Ho,{F,G)], 

where the random variables LLa[F,G) are recursively given by 

d 

H^,){F,G) = Y,5iG{^p%DF^), 

Ha{F,G) = H(^a^){F,H(^ai,...,a^:_i){F,G)). 

4.2. Conditional Malliavin calculus. In this section, we give the condi- 
tional version of some of the results established in Section 4.1. The proofs 
are very similar to the one-parameter case and are left to the reader. 

The first result is the conditional version of the duality relationship (4.1), 
which is the two-parameter version of [17], (2.12). 
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Proposition 4.4. Let s,t G with s <t. Let F he a random variable 
and let u he an adapted process such that E[/jg2 ds] < oo. Then 



a,2 



the following duality relationship holds: 



(4.4) E 



F 



Ur dWr 



[0,i]\[0,s] 



[o,t]\[o„ 



{DrF, Ur) dr 



The following norms are the two-parameter versions of those in [17], Def- 
inition 1. Let s,t G with s <t. For any function / G L^([0,t]",M'^), any 
random variable F G H)^'^ and any process u such that Ur G B'^'*' for all 
r G [0, t], we define 

= L2([o,f]\[o,s],K'^), 

1/2 



and 



IF 



\k,p,Hs 



{[0,t]\[0,s])' 



\\fir)\\i.dr 



n\F\f\J^,] + Y^E[\\D^F\f^^,\J^s] 



i/p 



u 



k,p 



r k ~| i/p 



Moreover, we write 7^ for the Malliavin covariance matrix with respect to 
Hs^t, that is, 

7^ = ((DFM)F^)h.,Ji<,,,<,. 

With this notation, we can state the following conditional version of in- 
equality (4.2), which is the two-parameter version of [17], (2.15). 

Proposition 4.5. Let s,t eM.\ with s <t. For any u G 'B^^'^'P{Hs,t), 
we have 

for some constant „ > 0. 



Finally, we give the conditional version of the two-parameter Burkholder 
inequality (see [18], A. 2, for the two-parameter nonconditional version). 

Proposition 4.6. Fixp > 1. There is a finite constant bp> such that 
for all adapted X = {Xt,t G M^) in L^{M.\ x $7) and all s,teM.\ with s <t, 



(4.6) E 



XrdWr 



[0,t]\lO,s] 



< bpE 



[0,t]\[0,s] 



X^dr 



p/2 



26 



R. C. DALANG AND E. NUALART 



4.3. Hyperbolic stochastic partial differential equations. Let 6, aj : — > M"^, 
^ ^ j ^ d, be measurable globally Lipschitz functions, where the vector- 
valued functions cti, . . . ,0"^ denote the columns of a matrix a = {o'j)i<ij<d- 

Consider the system of stochastic integral equations on the plane 

(4.7) Xi = XQ + y^[ a) {Xs) dWi + / h\Xs) ds, i G M^, 1 < i < d, 
~{Jm J[o,t] 

where the first integral is an Ito integral with respect to the Brownian sheet 
(as defined in [28], Chapter 4) and xq € M*^ is the constant value of the process 
Xt on the axes. It is well known (see [21], Lemma 3.1) that there exists a 
unique two-parameter, d-dimensional, continuous and adapted process X = 
{Xt,t G R^) that satisfies equation (4.7). In addition, E[sup^g[Q \Xj.\^] < oo 
for any p > 2 and t £ M^. In [21], Malliavin calculus is used to establish the 
following result. 

Theorem 4.7 ([21], Proposition 3.3). If the coefficients of a and b are 
infinitely differentiable with bounded partial derivatives of all orders, then 
Xi belongs to D°° for allteR\ andi = l,..., d. 

Assuming the latter infinite differentiability condition on the coefficients 
of a and 6, we state the following standard hypoellipticity hypothesis: 

Condition P. There is n > 1 such that the vector space spanned by 
the column vectors di, . . . ,cJd, cJjVcrj, 1 < i,j < d, cjj V((Tj Vcjfc), 1 < A; < 
d, . . . , cTj^ V(- • • (i7j„_i Vcjj„) • • •), 1 < ii, . . . , i„ < d, at the point xq is M'^, where the 
column vector aiVaj denotes the covariant derivative of aj in the direction 
of at. 

The following result uses Theorem 4.2 and gives the existence and smooth- 
ness of the density of Xt for any t away from the axes. 

Theorem 4.8 ([18], Theorem 2.4.2, and [21], Theorem 4.3). Under Con- 
dition V , for any point t away from the axes, the random vector Xt is nonde- 
generate and therefore has an absolutely continuous probability distribution 
with respect to Lebesgue measure on M!^. Moreover, its probability density 
function is infinitely differentiable. 

4.4. Gaussian-type upper bounds. In this section, we present some pre- 
liminary results and establish a Gaussian-type upper bound for the drift-free 
case with vanishing initial conditions. 
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Let X = {Xt,t G M^) be the unique solution of (4.7) with 5 = and xq = 0, 
that is, 

d „ 

(4.8) xi = Y aUX,)dWi, teM.l,l<i<d. 

We assume that the fohowing two hypotheses on the matrix a hold: 

Hypothesis PI. The coefficients of the matrix a are bounded and in- 
finitely differentiahle with bounded partial derivatives (we denote by T the 
uniform bound on the coefficients of a and its first partial derivatives). 

Hypothesis P2. Strong ellipticity.- ||o-(x)Cf = ELi(Ef=i > 
>0 for some p>0, for all x € M'^ and for all ^eR'^ with \\(\\ = 1. 

Note that Hypothesis P2 implies Condition P. Indeed, Hypothesis P2 
implies that the vector space spanned by the column vectors (Ti(x), . . . , a^ix) 
at any point x in M*^ is M'^, so Condition P holds. 

Fix s G M.^. Let Ps{io,-) be a regular version of the conditional distribution 
of the process {Xt — Xs,t G \ [0, s]) given J^g, as defined in Section 2. As 
in Theorem 4.8, we can check that under Condition P for P-almost all w G ^2 
and for any s <t, the law of Xt — Xs under Ps{uj,-) is absolutely continuous 
with respect to Lebesgue measure on W^. We let ps^t{^, x) denote the density 
oiXt — Xg under Ps{uj,-). We note that ps,t('^; 3;) is the conditional density of 
Xt — Xs given J^g. Therefore, for a random variable Y, we interpret E[y|.F5] 
as Es{lo,Y), where Es{uj, •) denotes the expectation under Ps{uj, •). However, 
for s G ]R+ fixed, there is Ns G J^s with F{Ns) = such that Ps,t(w,-) is 
defined ioi to £Cl\Ns and all s <t. 

Lemma 4.9. Assume Hypotheses PI and P2. Fix s,t G with s <t, 
and X G M'^. Let a be a subset of the set of indices {1, . . . ,d}. Then, for all 
Lv £^l\Ns and every a, 

Ps,t{uJ,x) = (-l)'^"''^'E[]l|xi_Xj>x%iGa, Xi-Xi<x\i(a}f^{l,...,d)(^t~^sA)\^s], 

where \a\ is the cardinality of a and, in agreement with Proposition 4.3, 
H(i^...4)iXt - Xs, 1) 

= 5(((7i;_xJ"'^(^t - Xs))% • • Sii^^l^ D{Xt - Xs)f) • • •)). 

Lemma 4.9 is a consequence of the integration by parts formula (4.3). 
It gives an explicit expression for the density of the process that will be 
very useful for further computations. When a = {1, . . . ,d}, the proof of 
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Lemma 4.9 follows along the same lines as the one-parameter noncondi- 
tional case (see [19], Corollary 3.2.1) and is therefore omitted. For the gen- 
eral case, see [22], (5.3), where a similar expression is obtained for the kernel 
of a stochastic semigroup. 

A key property is that the moments of the iterated derivatives of X are fi- 
nite; see [17], Lemma 6, where similar nonconditional estimates are obtained 
for one-parameter Brownian martingales in W^. 

Lemma 4.10. Assuming Hypothesis PI, for any < a < b < oo, p> 1 
and n>l, there exists a finite constant C > depending on a,b and the 
uniform bounds from Hypothesis PI such that for any s,t £ [a, 6]^ with s <t, 
on Q \ Ng, 

(4.9) 



sup 

zi<r,...,z„<r 
zi,...,z„,r£[0,t]\[0,s] 



E\\D 



Di':-\xir\^s]<c 



z\ 



for 1 < i,ki, . . . ,kn < d. 



Proof. We prove this lemma by induction on n. Suppose n = 1. For any 

z G [0,t] \ [0,s] and s,t € with s<t, the process {Di''\xi),l < k < d) 
satisfies the following system of stochastic differential equations (see [18], 
page 127): 

(4.10) Di'\Xl) = aliX,) + X: / Di'^a]{Xr) dW^ . 

■1 Jlz.t] 



Using Burkholder's inequality (4.6) for conditional expectations, for any 

Yn\Di>'\xir\:F,] 



i=l 



<Y2P-^\E[\ai{X^)nTs]+E 



J2f Df)a]{Xr)dW^. 



<2P-^ldTP + bpdP-^ J \Di^'^^j{Xr)\'^dr 



z,t] 
d 



p/2 



<2P~^\dTP + bpd^PbP~^TP f J2^M''^(^r)n^s]dr\. 

Finally, for z fixed, using a two-parameter version of Gronwall's lemma (see 
[18], Exercise 2.4.3) in the form 



(4.11) 



f{z,t)<A + Bf f{z,r)dr, 

J[z,t] 
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we conclude the proof of (4.9) for n = l. 

We now assume that (4.9) holds for n > 1. Apply n times the derivative 
operator to (4.10), which yields 

n+l 

1=1 



+ E 



ZlV---VZn + l,t] 



D^^-'^...Dllf\a]{Xr))dW}. 



For the first term, we use the induction hypothesis and the uniform bounds 
on the derivatives of the coefficients of a. For the second term, we use 
Burkholder's inequality and again the induction hypothesis and the bounds 
on the derivatives on the coefficients of a. We finally obtain 



i=l 



2lV---VZ„ + l,t] 



<Ci + C2E 

and the proof is completed using Gronwall's lemma. □ 



In the next lemma, we follow [17], Lemma 12, where similar estimates are 
carried out for one-parameter Brownian martingales in W^. 

Lemma 4.11. Assuming Hypotheses PI and P2, there exists a finite 
constant C > depending on a,b and the uniform bounds from Hypotheses 
PI and P2 such that, for any s, t E [a, 6]^ with s <t, on Q, \ Ng, 



(4.12) 



(E[{F(i,...,,)(Xi - Xs, 1)} Vs])'/' < C\\t - 



Proof. To simplify the notation, we write i/^^^ = H^i^ ^^-^{Xt — Xs,!). 
Using (4.5) and Holder's inequality for conditional Sobolev norms in Wiener 
space (see [29], Proposition 1.10, page 50), we obtain 



{l,...,d)\\0,2,Hs,t 



0,2,Hs 



< c\\H, 



t,s I 



'fU.Ell((7 



.s,t 
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Hence, to prove (4.12), it suffices to show that for each p>l and n = 1, . . . ,d, 
(4.13) \\D{xi - Xmi,,H., < cijt - s\\'/' 



for some finite constant p > and 



s,t 



l^ lit-, 11-1 



(4.14) Wiilx^-xj 
for some finite constant „ > 0. 



Indeed, (4.13) and (4.14) with n = 1 and p = 8 imply that 

II W*'* 11-^'= ^^^^1^211^ ^11-1/21117*.^ l|-^» 

llo,2,//,,t ^ cac c ||r s\\ ll^(i,„,,d-i) lli,4,H,,t- 
Iterating the process, we find 

\\^{l,...,d)\\o,2,Hs,t-^^"'^ ^ > ' 
which concludes the proof of (4.12). □ 

Proof of (4.13). Fix p > 1. By definition, 

\\D{xi - xDW-^"^^^^^^^ 

( n+1 

= n\D{xi-xi)r^jF,] + Y,n\DHxi-xi)\\i^,\j^s\ 



(4.15) 



i/p 



k=2 



Furthermore, for any z £ [0,t] \ [0,s], the process {Di^\xi - XI), l<k<d) 
satisfies the system of stochastic differential equations 

(4.16) Di'^^Xt - XI) = a\{X,) + f D^^^ai^X,) dW}. 

■ 1 J \z.t] 



By Burkholder's inequality (4.6) for conditional expectations and using the 
Cauchy-Schwarz inequality, for any f > 1, 



n\\D{xi-xi)r^^^\T. 



d 

E 

k=l 



\Df\xl-Xl)\Uz 



p/2 



'[0,t]\[0,^] 

< d^'/2-i(c||t - s||)P/2-i y E [ / \Df\x\ - X\)\^ dz 

U[0,t]\[0,s] 

<2P-idP/2-i(c||t_s||)P/2-i 
X J dTPc\\t-s\\ 
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<2P-i(c||t-s||f/2 

X I rfP/^TP + (c||t - s\\fl'^d^Pl'^-\{TY 
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dz 



X: sup E[|Z)«(X;)n^,]|. 



Finally, by Lemma 4.10, we get that there exists a positive finite constant 
ki{a,b,T) such that 



(4.17) 



E[\\D{Xl-Xl)rHj:Fs]<h\\t-s 



|P/2 



To estimate the second term on the right-hand side of (4.15), we apply n 
times the derivative operator to (4.16) to get 



n+l 

1=1 



+ E / Df^^)...Dt-^\^\a]{Xr))dW^. 

Proceeding as above, by Burkholder's inequality (4.6) for conditional ex- 
pectations, and using the Cauchy-Schwarz inequality and Lemma 4.10, we 
obtain, for all A; = 2, . . . , n -|- 1, 



(4.18) 



n\D'{xi-xi)r \T,]<k2\\t-sri\ 



where k2 is a finite positive constant. 

Finally, substituting (4.17) and (4.18) into (4.15) concludes the proof of 
(4.13). 

□ 

Proof of (4.14). Fix p > 1. By definition 



(4.19) 



k=i 



i/p 



32 



R. C. DALANG AND E. NUALART 



We use a standard argument to estimate the moments of the inverse of the 
Malhavin matrix. We follow [18], proof of (3.22), and [17], Lemma 10. Using 



the Cauchy-Schwarz inequality and Cramer's formula for (Txt-x^) ' 
can easily check that for all p > 1 , 

(4.20) '>^t-Xs 



IF, 



< c^Md^t^'xi-xX'^i^sV'^ X nmxt - x,)iix7^^i-F.] 

for some constant Cd,p > 0. For the second factor, we use (4.17) to get 

(4.21) E[||i?(Xi-X,)||^^ff ^Vs] <fci||*-s|P^^'"'^ 

for some finite constant ki{a,b,T) > 0. On the other hand, we write 



11/2 



4p(d-l)| ^ nl/2 



det 



s,t \ ■ r / T s.t 



-Xs 



IHI=iV^i^[M\[Ov 



i=l 



2 \d 

dz] . 



Using (4.16) and Hypothesis P2, for any h G (0, 1], we see that the expression 
in parentheses is bounded below by 

2 



E 



dz 



J2vJal{X,) + Y. I Df)a]{Xr)dW^. 
i=i \ j=i-'[^A 



> \Ap^ - h, 

where A denotes the area of the region [0, t — (1 — h){t — s)] \ [0, s] and 



d „ 



fc=l^[0,t-(l-h)(t-.)]\[0,s] 



X: / Df)a]{Xr.)dW}) 



dz. 



We choose y such that Ap^ = 4y ^/'^ and notice that since h <1, y > c:- 



t-s 



-2d 



In addition, as h varies in (0, 1], y varies in [c, oo) 



Applying Chebyshev's inequality for conditional probabilities, we find that 
for any q>2, 



< 



1 



1 



< 



1 



<F{h>y~'/'\J's}<y'^/''E[\h\'^\J^, 



Using Burkholder's inequality (4.6) for conditional expectations, for any 
g > 2, we have 

d 



iE[i4ri^s] < d^^"^M^^^"^^ E ^ 

i,j,k=l 



mms])- 



\Di^^a'j{Xr)\'^Urdz 
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By Lemma 4.10(i), the conditional expectation of the right-hand side is 
bounded above by some finite positive constant ^2(0, 6, T). Using the defi- 
nition of ^, we obtain 

Consequently, taking q> 2 + 2pd, 
E[(det7i:_xJ-'n-^.] 

= 2py2p-ip{(det7g„^J-i > y|^J dy 



< 



|t_5||2dp(5^)2dpp4dp 



roo 



00 



< + 2vk;^^^ I „2p-l-(g/d)+(2/d) , 

- lit - s\\^'ipibV2)^'ipp^'ip ^ y 

<h\\t-s\\-^'^p, 

where /C4 is a finite positive constant. Therefore, we have proved that 

(4.22) E[(det7j;_^J-2*'|.^,] < A;4||t - sH'^'^p. 
Substituting (4.21) and (4.22) into (4.20), we obtain 

(4.23) K[i((7i;-xj"').,ri-^.]<fc5(iit-sii-'^it-sf^('^-^))'/' 

= 1 1 ^ 5 1 1 ^ 

for some finite constant /C5 > not depending on i or j. This proves the 
desired estimate for the first term in (4.19). Turning to the second term, we 
claim that for all i,j = l,...,d and k = 1, . . . ,n, 

(4.24) mDHh'xl-xX'hfnf^l^^^ ^ h\\t - s\\-p 

for some finite constant fcg > 0. Indeed, by iterating the equality (see [18], 
Lemma 2.1.6) 

d 

k,l=l 

and using Holder's inequality for conditional expectations, we have 
supE[P^-((7l^!_xJ"').,ll>l-^.] 
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|P{r+l)|^il/(r+l) 



r=l fciH \-kr=k 

ki>l,l=l,...,r 



p(l+r)|^nl/(r+l) 



l^ |P(^+l)2|^nl/(r+l) 



where the supremum before the summation is over ii,ji,---,ir,jr G {1; • • • > c?}- 
By (4.23), for all i,j = l,...,d, 



E[|((7^'* 



■ I |p(r+l)2 
Xt-Xj Ail 



~p{r+l)2 



for some finite constant > 0. 

For the other factors, express D'^ilxl-Xs^^j using the definition of 'Jxl-Xs 
and use the Cauchy-Schwarz inequality twice and (4.18) to get 



[o,t]\[o. 



Dr{Xl-Xl).DriXl-X^,)dr 





p 




) 











1=0 



E 



[o,t]\[o„ 



D^Dr{Xl-Xi 



D''-^Dr{Xl -Xi)dr 



p 








TT®k 

"-a,t 



<kdt-s\\P 



a/2 



x{E[\\D'-^D{X}-Xi)\r^,,_,^,,\Ts]r^} 



for some finite constant /cg > not depending on i and j. This concludes 
the proof of (4.24). Finally, substituting (4.23) and (4.24) into (4.19), we 
conclude the proof of (4.14). 

The proof of Lemma 4.11 is now complete. □ 



The following result is a consequence of Lemmas 4.9 and 4.11. 



Lemma 4.12. Assuming Hypotheses PI and P2, for any < a < b < oo, 
the density function pxt {x) is uniformly hounded for t G [a, and x . 
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Proof. Using Lemma 4.9 with s = and a = {1, . . . ,d} and the Cauchy- 
Schwarz inequahty, we obtain 

pxA^)<{E[{H^,_,^{Xt,l)}Y^'. 

By Lemma 4.11 with s = 0, there exists a finite positive constant C2 de- 
pending on a, b and the uniform bounds from Hypotheses PI and P2 such 
that 

E[{i/(i,...,d)(Xi,l)}2]<C2, 
which proves the lemma. □ 

The next proposition is the main result of this section. 

Proposition 4.13. Assuming Hypotheses PI and P2, for any < a < 
b < cxD, there exists a finite positive constant c depending on a,b and the 
uniform bounds from Hypotheses PI and P2 such that for any x G M.'^, s,t & 
[a, b]'^ with s <t and for co \ Ng, 

(4.25) PsA^^x) < c\\t - s\\-'^/'^ exp( I^^V 

V c||i — 

Proof. Apply the Cauchy-Schwarz and Holder inequalities for condi- 
tional expectations to the expression of Lemma 4.9 with a = {i £ {1, . . . , d} : > 
0} to find that 

(4.26) p.,t(^,x) < f[{n\Xt-Xi\ > \x'\ \Ts})'/(''\E[{Hll^^^^^f]f'. 

i=l 

Consider the one-parameter martingale [M„ = (M^, . . . , M^),0 < u < |t| — 
\s\] defined by 



d 



[ f^j(^r-) dW^, if < u < ti - Si, 

j^lJ[si,si+u]x[0,S2] 

K-.i+E / a){Xr)dW^, 

if ti — si < u < \t\ — \s\, 

for alH = 1, . . . , d, with respect to the filtration (^u, < u < |t| — |s|) defined 
by 

G = /-^(^i+^.^a)' ifO<u<ii-si, 
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Notice that Mq = 0, M\ 



Xt-Xs and go = J^s. By [1], (2.9) 



d „ 

= E/ ^j(^r) 



Moreover, Hypothesis PI and the Cauchy-Schwarz inequaUty imply that 
(M*)|t|„|^[ < C\\t - s\\, where C = bd2^/^T^ for all I < i < d. Applying the 
exponential martingale inequality [18], A. 2, we get 



Finally, substituting (4.27) into (4.26) and using Lemma 4.11 concludes 
the proof of (4.25). □ 

4.5. Gaussian-type lower bounds. In this section we present a lower bound 
of Gaussian type for the density of the random variable Xf for any t away 
from the axes, where X = {Xt,t S M^) denotes the solution of (4.8), and 
we present an analogous lower bound for the conditional density of Xt — Xs 
given J^s, when s <t. These results are an application of results by Kohatsu- 
Higa [13], where Gaussian-type lower bounds are obtained for the density 
of a general class of uniformly elliptic random variables on a Wiener space, 
which generalize the lower bound estimates for uniformly elliptic diffusion 
processes obtained by Kusuoka and Stroock [14] . 

Theorem 4.14. Assume Hypotheses PI and P2 and let X = {Xt,t £ 
M^) be the solution of (4.8). Then for any < a < b < oo, there is a constant 
C which depends only on a,b and the uniform bounds from Hypotheses PI 
and P2 such that, for all s = (si,S2) £ [oi&]^; 



Proof. Following [13], Theorem 5, we need to prove that for fixed s = 
(si,S2) G [0)^]^5 the d-dimensional random vector Xs is uniformly elliptic, 
with constants that do not depend on s. Following the same notation as 
in the Main setup of [13], we replace their [0,r] x A with [0, si] x [0, S2] 
and we set g{-,-) = 1. The underlying one-parameter filtration is defined as 
{Tl,0<t< si), where 



(4.27) F{\Xl-Xl\>\x'\\Ts}<2exp 



2Cllt- s 



X' 



) 



l<i<d. 




J^l = (7{W{ZI, Z2),{ZI, Z2) £ [0,t] X [0,S2]}. 



(In this proof, t denotes a real number, as in [13].) 
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Consider a sufficiently fine partition {0 = Iq < • • • < tj\f = si} of the in- 
terval [0, si]. To simplify the notation, we write Rt„_i,t„ for the rectangle 
[tn-i,tn] X [0,S2]. For any 6 = {01,62) G Rt„-i,t„, we write irO for the point 
^2)) that is, the orthogonal projection of 9 on the vertical line through 
(t„_i,0). We write i?7re,6» for the rectangle [tn~i,Oi] x [0,02]- Finally, for any 
9 £ Rtn-i,tn and r] S RnOfi, we write ^ A 77 for the point (771,^2)) that is, the 
orthogonal projection of rj on the horizontal line through Tr9 and 9. Let 
An-i{g) be the quantities denoted Aj_i(g') in the Main setup of [13], that 
is, 

An-l(5)= / \\g{t,-)\\^2(^[Q^,^]^^d)dt={tn-tn-l)s2 

for 71 = 1,..., iV. Note that ||c/|Il2([o_5^]x[o,s2],iR'*) = '^i^a- 

We now define the sequence of nondegenerate random vectors Fn required 
in the Main setup of [13] by F„ = Xf^^g^^ for < n < A^, with = Xsj^^s2- 
Notice that F„ = , . . . , and F^ = . By (4.8) , 

(4.28) Fi - = Y^\ a]{Xe) dW^ l<i<d. 

The first objective is to find the Ito expansion of order / > 1 of the ran- 
dom variable Fl^ — F^_i to obtain the approximations F^ required in [13], 
Theorem 5. For this, we need to introduce some notation. We define the 
multiindex /? S Un>i{0) with length l{f5), and write f3 = (/3i, . . . ,/3;(^)). 
We write — /? and (3— for the multiindex obtained by deleting the first and 
the last component, respectively, of the multiindex /?. For completeness, we 
write {v} for the multiindex of length zero. These multiindices will be used 
to write multiple integrals, some of which are stochastic and some are de- 
terministic: when the index is 1, d^Wg denotes dWg , and when the index 
is 0, d^Wg denotes d9. Usually, these integrals are not taken on the whole 
space [0,si] x [0,^2], but on subsets of it. The subset most often is of the 
type [t„_i,t„] X [0,52]- In general, these integrals are denoted by //3(/i/3) for 
an JP"/^^ -measurable random process h/^ such that 

sup ll^/3|lL2(([J„_i,t„]x[0,S2])'('3),Krf)(^) < 

Given a family of M functions fm ■^'^ ^ M, 1 < m < M, and M points 
9^,..., 9^^ placed on the same vertical line in i?4„„i,t„, we define for any 
1 < r < d and ^ € Rtn-i,t„ the following operations on random variables: 

\m=l / 
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M / M 



M / M \ d Qj: 

m=l \n=l,n^m /k=l ^ 



M 



V m=l 



\m=l \n=l,n=^m / 



d 



_ d'f n 
k,l=l ' 



+ E 

(0 n fn{Xe«Ai)] 

Lmi,m2=l \n=l,nj^mi,m2 / 

Note that the points 6^ A ^, . . . ,0*^ A ^ and ^ are also placed on a single 
vertical line in Rt„_i,tn- 

These operations will be used to apply the standard multidimensional Ito 
formula to f{Xgi , . . . , Xqm), when f{x^, . . . , x^^) is of the form 11™=! 
where the fm are as above, 1 <m< M. Indeed, 

""1 ^ (X^^ Qi,. . .,X^_^ qm), tn-l <Ui< tn, 

is an Md-dimensional martingale with mutual covariation 
Therefore, according to the Ito formula [2], 

M M 
Y\_ fm{Xgm) = Y[ fm{X-^em) 
m=l m=l 



d „ I M 

(4-29) +E / n /-(^e^ 

r=l •^■^tn-l.in \m=l 

d „ / M \ 

+ El n/-(^^-) Me. 

r=l-^-^*n-i^tn \m=l / 

We are now going to iteratively apply the Ito formula to the integrands in 
(4.29), to obtain an Ito expansion similar to the one presented in [12], Chap- 
ter 5. For this, we introduce some additional notation. Given a multiindex 
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j3, we define the functions 

%(0) = c7j(X,), if/J = W, 

%,,(e,0 = iW(cTi(Xe)), if/3 = (l), 
Wp^^^,{e,i) = L%^^{a]{Xe)). if/3 = (0), 

wliere 

We define the multiple (Ito stochastic/deterministic) integrals by 



= E 

(4.30) i,n,-,n(^)-i=i 



.,n(/3) 



where the domain of integration of each integral is given by the indicator 
functions that appear in the definition of the operators LP and L^. 
Note that the functions j, ^ , . . . ,r, _ 1 ' ^ ^ • • • > ^' ^'^^ ~ ^ ) 

are sums of func- 
tions of the form Jlri^i fniXgn ), where rM*^ — > M are derivatives of some 
order of the coefficients of the matrix a, and 6^ , . . . , 9^^^'^ are l{f3) points 
placed on a single vertical line in Rt„_^,t„ that depends on 9,^,^, . . . . 
We now define multiple Ito stochastic integrals I^, by modifying the in- 
tegrand in the definition of each point in (4.30) is projected onto the 
vertical line through {tn-i,0), that is, 

= E i (/, (■^■(/En/»(^-) 



X d^m-iwp ] dWi 
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Lemma 4.15. The ltd expansion of order I > 1, 

holds, where the sets Ai and Bi are recursively defined by Ai = {(1)}, Si = 
{(1,1), (0,1)}, Ai+i = {Ai,Bi} andBi+i = {(3:-(3eBi} forl>l. 

Proof. We prove this lemma by induction on /. By the standard multi- 
dimensional Ito formula, applied to the random variable aj{Xg) with respect 
to the first coordinate, it follows from (4.28) that 



= E/ <^j{X^9)dWi 



for all l<i<d, which proves (4.31) for / = 1. 

Now assume that (4.31) holds with Ai-i and Bi^i for / > 1. Then we apply 

the Ito formula (4.29) to all the random variables of the form nn=i fn{Xgn) 
that appear in the stochastic integrals of the second term of (4.31) with Bi-i. 
Here, M = /(/3), the functions — > M, 1 < m < /(/3), are derivatives of 
order greater than or equal to of the coefficients of the matrix a and 
9^,. . . are Z(/3) points placed on a single vertical line in Rt„_i.tn - The 

first term on the right-hand side of (4.29) gives a new IJ^ integral and the 
two other terms each give a new I/j integral. Note that integrals with the 
operator add a 1 to all the multiindexes of Bi-i and that integrals with 
the operator add a 0. This yields (4.31) for I with Ai = {Ai^i,Bi^i} and 
Bi = {P:-PeBi_i}. □ 

Continuing the proof of Theorem 4.14, with this result we define, follow- 
ing [13], Theorem 5, the approximation of order / > 1 by 

F; = ((t„-t„_i).2)^'+'^/'^;+F^i 

+ E^^^[^^-.,n,...,.(«_.(^>e\...,e'(^)-^)]^.„ 
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where n = 1, . . . , A^) is an i.i.d. sequence of d-dimensional N{0,I) ran- 
dom variables independent of the Wiener process W. In the setting of [13], 
Theorem 5, we have k = l and set 



To simphfy notation we write /i^ for 

^/3-,j,ri,...,n{^)_i(6',?\---,C'' 

With these definitions and since by Hypothesis PI, for any multiindex /3 G 
A, 

sup ll^/3llL2(([i„_i,t„]x[0,S2])'('3),Kd)(^) < 

(HI) of [13], Theorem 5 is satisfied. To prove (H2a) of [13], Theorem 5, note 
that by definition 



'"TL 

n ^'n\\n,p 



and, for each f3 £ Bi, 



n,p 



(4.32) 



= |iE[i//3[/^j5]i?.„_,,„ri-^Lj 



+ 



1/p 



fit. 



IIP 



where Ht„_^,tn = L'^{[tn-i,tn] x [0,S2], 

Fix P (z Bi and recall that /(/3) =1 + 1. To estimate the first term in (4.32), 
we use Burkholder's inequality (4.6) for conditional expectations and the 
uniform bounds on the derivatives of the coefficients of a to get 

(4.33) E[|/^[/.^]«,_^,,„ f\rl_J < C{{tn - tn-l)s2f^'^''^ 

for some positive finite constant C > independent of the partition, ti, s 
and cj G il. 

For the second term in (4.32) fix j E {l,...,g}. Using formula (1-46) 
from [18], we can easily check that the term Z)-^/^[/i^]ij(^ contains mul- 
tiple (Ito stochastic/deterministic) integrals of orders / + 1,...,/ + 1— j. 
Therefore, again using Burkholder's inequality for conditional expectations 
and the uniform bounds on the derivatives of the coefficients of o", we obtain 



(4.34) E 
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for some positive finite constant C > independent of the partition, ti, s 
and Lo gQ. 

Finally, using (4.33) and (4.34), we get that there exists a constant C 
independent of the partition, ti, s and w G Jl, such that 

11^^ - K\\n!p-' < C{{tn - tn-l)s2f^^^'\ 

which proves (H2a) of [13], Theorem 5, with 7 = 1/2. 

Now, using exactly the same argument that led to (4.22) we can easily 
check that 

(E[(det7j.:--*")-n.^,Lj)'/^ < C((i„ - t„_i).2)"', 

where 7f„"''*" = ((^^n> ^^/i)L2([t„_i,t„]x[o,s2],K'*))i<ij<rf> for some positive 
finite constant C independent of the partition, ti and uj ^Vt. Therefore, 
(H2b) of [13], Theorem 5, is proved. 

Condition (H2c) of [13], Theorem 5, can be rewritten in our case as 

2 

C2<{{tn-tn-l)s2) ^ ' ~'~ ■ 



for ^ G M'^ with ||^|| = 1. This is obviously satisfied by Hypotheses PI and 
P2. 

Finally, condition (H2d) of [13], Theorem 5, is satisfied using the same 
argument that led to condition (H2a), because the higher-order integrals are, 
in the J~t„_i conditional D"'^ norm, smaller than (tn — tn-i)s2- Theorem 4.14 
is now proved. □ 

The next result is the conditional version of Theorem 4.14. Let ps,t{'^-,x) 
be the density of Xt — Xs under the conditional distribution Ps{ui, •) defined 
just above Lemma 4.9, for w G $7 \ Ng. 

Theorem 4.16. Assuming Hypotheses PI and P2, for any < a <b < 
00, there exists a positive finite constant c, which depends on a,b and the 
uniform bounds from Hypotheses PI and P2, with the following property: 
For all s, t G [a, 6]^ with s <t, x G M"^, and for u) e^\Ns, 

Ps,t{^,x) > c\\t - s||"'^/^exp (- |[^^^ ) . 

V cjji — 

Proof. Fix s,t G [a, 6]^ with s <t. As in Theorem 4.14, Theorem 4.16 
is an application of [13], Theorem 5. We replace their [0,T] by the union of 
two segments with extremities (si,S2), (^1,^2) and (ti,S2), (ii,i2)) and we 
set (?(•,•) = 1. We consider a sufficiently fine partition {s = to < ■ ■ ■ < tjy = t} 
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of the union of the two segments, where ti = {tj,tf), 1 < i < N . Note that 
there exist two constants ci{a,b) and 02(0,6) such that 

Ci ^\\tn — tn~l\\ < A„_i((7) < Ci||t„ — tn-l|| 

and 

C2 - s|| < ll5llL2([0,t]\[0,s],R<*) < C2||t - s\\. 

In this case, ah of the Wiener stochastic integrals are taken on sub- 
sets of one of the [0,tn\ \ [0,i„_i]. Then, the projection of a given point 
9 E [0,tn] \ [0,tn-i] depends on its position. Namely, if 6 is in the rectan- 
gle [si,ti] X [0,^2], 7^0 will be its orthogonal projection on the vertical line 
through (t^_]^,0), while if 6 is in the rectangle [0, ti] x [525^2]) will be its 
orthogonal projection on the horizontal line through {0,t'^^_i). 

Define = - for 1 < i < d, < n < N. Then 

d „ 

Fl-Fl_^ = Y, a){Xe)dWi, l<i<d. 

r;^i[o,i„]\[o,t„-ii 



J- 



Now, using the same notation as in the proof of Theorem 4.14 and pro- 
ceeding exactly along the same lines, we can easily establish the following 
result: 



Lemma 4.17. The ltd expansion of order I > 1, 

Fl - F;_i = ^ /^[/l^j.^j. . . • ,C«(/3)-l)][0,t„]\[0,t„-i] 

+ Yl ^^[^/3-J>i,...,r;(/3)_i(^'^l'---'^H/3)-l)][0,t„]\[0,t„_i]' 

holds, where the sets Ai and Bi are recursively defined by Ai = {(1)}, Bi = 
{(1,1), (0,1)}, Ai+i = {Ai,Bi} andBi+i = {(3:-(3GBi} forl>l. 

Again, following [13], Theorem 5, we define the approximation F^ of order 
/ > 1 by 

— llf -f J|('+l)/2 7« I Z?* 
n ~ \\'"n ''n— 111 ^n'-'-n—l 

+ E ^^[^/^-J,n,...,n(/3)-l(^'^l' - ■• '^K/3)-l)][0,i„]\[0,i„-i]' 

where (Z^, n = 0, . . . , A^) is an i.i.d. sequence of d-dimensional N{0,I) ran- 
dom variables independent of the Wiener process W. 

The remainder of the proof follows exactly as in the proof of Theorem 4.14. 
We note that the constant c in the conclusion of Theorem 4.16 is also uniform 
in uj £Q\Ns, s,t£ [a, 6]^ with s<t and x G M''. □ 
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5. Potential theory for hyperbolic SPDEs. In this section, we extend the 
results obtained in Section 2 to the solution of equation (4.7). The proofs 
make use of Malliavin calculus and are an application of the results of Section 
4, which contains the technical work. 

5.1. The case b = 0. Let X = {Xt,t£ R\) he the unique d-dimensional 
adapted continuous process defined on {Q,Q,F) that solves (4.8). The aim 
of this section is to establish the following result. 

Theorem 5.1. Assuming Hypotheses PI and P2, for all < a <b < oo 
and M > 0, there exists a positive finite constant K depending on a,b,M,p 
and the uniform bounds on the coefficients of a and its derivatives, such that 
for all compact sets Ac {x ^M.'^ : \\x\\ < M}, 

K-^ Caprf_4(^) <F{3t£[a,bf:XteA} <KCapa^^{A), 

where Cap^ denotes the capacity with respect to the Newtonian (3 kernel 
kf3{-) defined in Theorem 3.1. 

Proof. To prove Theorem 5.1 it suffices to prove that Hypotheses Hl- 
H3 of Theorem 2.4 hold for the process X. Since Hypothesis H3 is only 
used for the upper bound in the statement of Theorem 5.1, we only prove 
Hypothesis H3 when (i > 4, since the upper bound is trivially satisfied for 
d < 3. Fix < a < 5 < oo, M > and recall that Hypotheses PI and P2 are 
satisfied. 

Verification of Hypothesis HI. Fix x G M"' such that < M. 
Using Theorem 4.14, we find that 

/ PxAx)dt>C I (M2)~''/'expf-Mi)di 

>C(6-a)V^exp(^-^), 

which shows that Hypothesis HI is satisfied. 

Verification of Hypothesis H2. Fix x and y such that \\x\\ < M 
and llyll < M. Let s,t G and assume that s <t. Clearly, 

PXt,Xs {x, y) = PXt-Xs,Xs {x - y, y) 

= Pxt-Xs\Xs=yix - y)pxs{y)- 

Let z = X — y and let E he a Borel subset of M'^. By Proposition 4.13, we 
have, a.s., 

F{Xt - Xs e E\Ts} < I c\\t - s\\-'^/'^ e-x.J — \ dz 

JE \ c\\t — s\\J 
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for some finite constant c > 0. Take tlie conditional expectation of both sides 
with respect to cr{Xs) to find that 

Px,~x^\x.,=y{z) < c\\t - s||-'^/2exp(^- J^^) . 

By Lemma 4.12, PXsiu) is uniformly bounded over s G [a,b]'^ and x G W^. 
Therefore, for all s,t G [a, with s < t, we have proved that 

c\\t — s\\ 



(5.1) px„xAx,y)<c'\\t-s\\-''^'exp 

for some finite constant c' > 0. 

We now assume that s,t £ [a, 6]^ with s and t ^ s. Let Ei and E2 be 
two Borel subsets of W^. Using the conditional independence property, we 
obtain 

F{Xt £ Ei.Xs £ E2} = E[F{Xt £ Ei,X, £ ^zl^sAt}] 

= E[F{Xt £ E^lJ'.^.tmXs G E^lJ'sAt}] 

= E 



El JE2 

Set p{c,x) = c"'^/^ exp(— ||xp/(2c)). By Proposition 4.13 and Lemma 4.12, 
the right-hand side is bounded above by 

CE / dxp{c\\t-{sM)\\,x-XsM) I dyp{c\\s-{sM)\\,y-XsM) 

JE2 



= C dx dy dzp{c\\t — {s At)\\,x — z)p{c\\s — {s At)\\,y — z). 

JEi Je2 JR'i 

For x fixed, do the change of variables u = z — x and use the fact that 
p(cj^,-) is an even function to write the inner integral as a convolution of 
two Gaussian densities, which is equal to 

p{c\\t - {s At)\\ + c\\s - {s At)\\,y - x) < Cp{c\\t - s\\,y -x) 

by the triangle inequality. We conclude that 

F{Xt £Ei,Xs£E2}<C [ dx [ dyp{c\\t 

JEi JEo 

and, therefore, 

(5.2) (^,y)<C7||t-s||-^/2expf-^^''^ 

V c\\t — s 

This shows that an estimate like (5.1) holds also when s and t ^ s. 



s\\,y-x) 

lEi JE2 
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Using (5.1) and (5.2), we get 

/ / PXt,Xs{x,y)dtds 

d/2 ^.,J Ik -2/ "^ 



<C / \\t- s\\-''l'e^Y>[-\ — "-\]dtds. 

J[a,fe]2 J[a,fe]2 V C.\\t-s\\J 

Fix t and use the change of variables u = t — s to see that this expression is 
less than or equal to 

4C{b-af f \\u\r''/\xp( -^^^^-^] du. 

J[0,fe-a]2 V C\\u\\ ) 

Next, use the change of variables m=||x — y|pc~^2;, to see that this is less 
than or equal to 

c'^/2-24c(5_a)2||2._y||-'^+4 /■ ||z||-'^/2exp(-l/||z||)(iz. 

J[0,c(6-a)/||a;~y||2]2 

When d > 4, applying Lemma 3.5 with /? = djl > 2 and a = 1/2 in the same 
way as in the proof of Theorem 3.1 shows that Hypothesis H2 is verified 
with A;(-) = A:rf_4(-) and d > 4. When d < 4, the above expression is bounded 
and Hypothesis H2 holds with k{x) = \. 

Verification of Hypothesis H3. Assume d > 4. Fix x G M*^ with 
\x\ < M. Use the Gaussian-type lower bound for ps^t{^,x), obtained in 
Theorem 4.16, to see that for all s G [a, 6]^ and for almost all cu, 

Ps t{^,x) dt> c [ \\t — s\\~'^^'^ expf — J — ^- — dt 

[fe,26-a]2 ' J[6,2b-a]2 V c\\t - s\\ J 

for some finite constant c > 0. Now, using the change of variables z = c{t — s)/ 
we see that this expression is 

>c^/2-2||^||-d+4 f ||z||-'^/2exp(-l/||z||)(iz. 

J[0,c(6-a)/||x||2]2 

Finally, Lemma 3.5 with (3 = d/2 > 2 and a = 1/2 shows that Hypothesis 
H3 holds for d > 4 in the same way as in the proof of Theorem 3.1. 

The proof of Theorem 5.1 is complete. □ 



/ 

J\b. 



5.2. The case 6^0. The aim of this section is to extend Theorem 5.1 to 
the case 6^0. Our main tool is Girsanov's theorem for an adapted transla- 
tion of the Brownian sheet (see [20], Proposition 1.6). 
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Let Y = (Yj, t G R^) denote the d-dimensional adapted continuous process 
defined on {0,,Q,¥) satisfying (4.7) with xq = 0, that is, 

(5.3) Y; = y [ aUYs) dWi + / h\Ys) ds, teRl,l<i< d. 

We introduce the following condition on the vector b: 

Hypothesis P3. For some constant N, for all 1 < i < d and x G M*^, 
\¥{x)\ <N. 

Consider the random variable 

Lt = exp\- [ a-\Ys)b{Ys) ■ dWs - I f \\a-\Ys)b{Ys)\\l, ds 
I J[o,t] J[o,t] 

We have the following Girsanov theorem. 

Theorem 5.2 ([20], Proposition 1.6). The random variable Lt is such 
that E[Lt] = 1. If¥ denotes the probability measure on {fl,J-') defined by 

then Wt = Wt + /[q j] a~^{Ys)b{Ys) ds is a standard Brownian sheet under P. 

Consequently, the law of Y under P coincides with the law of X under P, 
where X = {Xt,te 1R+) is the solution of (5.3) with 6 = 0. 

The following result is the extension of Theorem 5.1 to the case 6^0. It 
is sufficient to characterize polar sets of Y . 

Corollary 5.3. Assuming Hypotheses P1-P3, for all < a <b < oo 
and e,M > 0, there exists a finite positive constant depending on a, b, e, 
A, M ,N , p and the uniform bounds of the coefficients of a and its derivatives, 
such that for all compact sets A<Z {x ^ • ll^ll < -^}; 

A7nCap,_4(A))^+^ < P{3t E [a, bf ■.Yt(^A}< K,{C3.Vd-M)f'^^" ■ 
Proof. Consider the random variable 



Jt = exp 



/ a-\Xs)b{Xs)-dWs + \ I \\a-\Xs)b{Xs)\\l,ds 



Fix < a < 6 < OO and let Gx = {3 1 G [a, bf:Xte A] and Gy = {3 1 G 
[a, : Yt £ A}. By Theorem 5.2, 

(5.4) F[GY]=MlGy]=^p[tGyLt^]=mtG^Jt~']- 
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Let e > and apply Holder's inequality: 

Rewriting the last inequality we obtain 

Let r > 0. By the Cauchy-Schwarz inequality 



exp 



X E, 



/ {-2)ra-\X,)b{Xs)-dWs 
.J[o,t] 

-\ f Ar^\\a~\XMXs)\\l,a 

[ i2r' + r)\\a-\XMXs)\\l.ds 
J\o.t] 



1/2 



The first expectation on the right-hand side equals 1 since it is the expecta- 
tion of an exponential martingale with bounded quadratic variation (see [9], 
Chapter 3, Proposition 5.12). By Hypotheses P2 and P3, the second factor 
is bounded by some positive finite constant. Therefore, 

Ep[j//"] < ke 

for some constant > 0. Finally, by Theorem 5.1 there exists a positive 
finite positive constant K such that P[Gx] > Cap^„4(A), which concludes 
the proof of the lower bound. 

The upper bound is proved along the same lines. Let e > and apply 
Holder's inequality to the right-hand side of (5.4): 

P[Gy] =Ep[lGx'/r'] < (P[Gx])'/'+^(Ep[jr^'+'/'^])^/'+^ 
Let r > 0. Again by the Cauchy-Schwarz inequality we find that 



1 < (ep 


exp 


/ 2ra 






-Jlo,t] 



[0,t] 



Ar^\\a-\Xs)b{Xs)\\l, ds 



1/2 



X E 



exp 



[o,t] 



i2T''-r)\\a-\XMXs)\\l.ds 



1/2 



The first expectation on the right-hand side equals 1 since it is the expec- 
tation of an exponential martingale with bounded quadratic variation, as 
above. By Hypotheses P2 and P3, the second factor is bounded by some 
positive finite constant. Therefore, 

Ep[jr^'+^/^^] < K 
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for some constant > 0. Finally, by Theorem 5.1, there exists a positive 
finite positive constant K such that P[Gx] < -f^Cap^_4(^), which completes 
the proof of the Corollary 5.3. □ 

As a consequence of Corollaries 2.5 and 5.3 we obtain the following ana- 
lytic criterion for polarity for the process Y: 

Corollary 5.4. Assume Hypotheses P1-P3. Let E he compact subset 
of . Then E is a polar set for Y if and only if Cap^_4(£') = 0. 

Finally, Theorem 3.3 and Corollary 5.3 give the stochastic codimension 
and Hausdorff dimension of the range of the process Y: 

Corollary 5.5. Assuming Hypotheses P1-P3, 
codim{y((0, +oof)] = {d- 4)+ 

and if d> A, then 

dim-H{y((0,+cx))2)} = 4 a.s. 

5.3. Critical dimension for hitting points. For an A^-parameter R'^-valued 
Brownian sheet W = {Wt,t G M:!J^), Orey and Pruitt [23] showed that for any 



1, \id<2N, 
0, ifd>2iV. 



This fact also can be obtained as a consequence of work by Khoshnevisan 
and Shi [11]. Corollary 5.4 immediately yields an analogous result for the 
solution Y of (5.3): 

Corollary 5.6. Assuming Hypotheses P1-P3, for any x G , 
V{3t^Mi^:Yt = x}>{) if and only if d < A. 
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